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The direct detection of dark matter has still not been achieved. Many attempts are 
underway including the global DarkSide-20k collaboration. The light-sensing technol- 
ogy (SiPMs) used in this experiment will be tested at Warwick. This paper documents 
the building of model circuits in Multisim which constitute a toolkit for investigation 
of this technology. The typical SiPM output was simulated in order to prepare for 
the testing process. The simulation was realistic and could be effectively analysed. 


The results of this work will prove useful for the testing process, the realisation of the 


DarkSide-20k experiment, and will provide tools for wider research. 


I. Introduction 


One of the latest endeavours in the direct 
detection of dark matter is DarkSide-20k. 
So-called as it will use 20 thousand kilo- 
grams of detector material. This global col- 
laboration aims to detect WIMP dark mat- 
ter through the use of a liquid Argon detec- 
tor. The detector relies on the production 
of a light signal. The light sensors used in 
this detector will be Silicon photomultipli- 
ers (SiPMs). SiPMs have not been used in 
a detector such as this before. They offer 


numerous benefits over the current technol- 
ogy - photomultiplier tubes (PMTs). 

This document details a project which had 
the aim of developing a computer-based 
model of a SiPM and using this to cre- 
ate analysis software. This software will 
be used to test and characterise the SiPMs 
before they go into the DarkSide-20k detec- 
tor. This process occurred in three stages. 
Firstly, a series of model circuits were built 
which output a SiPM signal. Secondly, a 


Monte Carlo method was used to simulate 


different types of noise that occur in SiPM 


signals before finally writing software to 
analyse it. 

This paper will outline the concepts behind 
dark matter detection and SiPMs before 
moving on to describe the work that was 


undertaken. 


A. Evidence For Dark Matter 


Despite its existence first being proposed 
in 1932, dark matter didn’t come into the 
fore until the 1980s. Jan Oort (1932) and 
Fritz Zwicky (1933) were the first to spot 
that galaxies seem to have a stronger grav- 
itational pull than their luminous masses 
would suggest [1, 2]. Since then, studies of 
the rotation velocities of galaxies, gravita- 
tional lensing and other phenomena have 
provided a plethora of indirect evidence for 
some undetected mass in the universe - 
dark matter [3, 4, 5, 6, 7]. It is estimated 
that 83.8% of all matter is dark [8], yet di- 
rect detection remains elusive. The relative 
proportions of the different mass and en- 
ergy components of the universe are shown 
in Figure 1. 

Dark matter has (thus far) only been ob- 
served through its gravitational effects. 
Therefore, it follows naturally that dark 
matter cannot (or can only very rarely) in- 
teract with normal baryonic matter. The 
question of what this invisible mass is com- 
posed of remains one of the greatest chal- 


lenges of modern physics. Many theo- 
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Figure 1: The relative abundances of the differ- 
ent constituents of the universe. Over 95% of the 
universe is made up of dark matter and dark en- 


ergy [9]. 


ries have been proposed, but there is still 
no solution to the problem. Dark mat- 
ter cannot be directly studied yet so, to 
build a testable model, grand assumptions 
must be made about its composition and 
properties. A common assumption is that 
dark matter consists of particles that are 
yet to be directly detected. Some of the 
particle-nature candidates for dark matter 
are shown in Figure 2. 

An early explanation for the observed grav- 
itational discrepancies were dim stars, stel- 
lar remnants and primordial black holes 
(PBHs) [10]. Studies have since shown that 
these objects could only make up a minute 


proportion of dark matter [11, 12, 13]. 


B. The WIMP Model 


It is now known that dark matter consists 
of non-baryonic material. The two most 
popular candidates for dark matter parti- 
cles are axions and weakly interacting mas- 
sive particles (WIMPs). Both of these had 


already been proposed in particle physics, 
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Figure 2: The different dark matter candidates and their predicted mass ranges. WIMPs (the focus of 


the DarkSide-20k project) are highlighted in blue [8]. 


independent from the search for dark mat- 
ter. The focus of this paper is on WIMPs 
so axions will not be discussed here. The 
motivations for and theory behind axions 
can be found in [14]. 

WIMPs have been proposed in theories go- 
ing beyond the standard model. They are 
defined as massive, electrically neutral par- 
ticles which only interact through the weak 
force. It is assumed that WIMPs were 
thermally produced in the early universe 
and annihilate with each other. Taking 
into account that they must have cross- 
sections on the scale of the weak force, the 
number of WIMPs remaining in the uni- 
verse today can be estimated. The result 
states that WIMPs should make up approx- 
imately 25% of all matter and energy in the 
universe - just below the observed abun- 
dance of dark matter (see Figure 1). This 
coincidence is known as the “WIMP mir- 
acle” and makes this particle one of the 


frontrunners in the search for dark mat- 


ter [8, 10]. 


C. WIMP Detection and Detector 
Types 


Due to their non-interactive nature, 
WIMPs have not yet been detected. There 
are four ways in which the search for 


WIMPs is proceeding: 
e WIMP production at CERN. 


e Direct detection from scattering exper- 


iments. 


e Indirect detection from dark matter 


annihilation products. 
e The search for dark stars. 


See [10] for descriptions of the particle de- 
tection attempts and see [15] for a review of 
dark stars. This paper will focus on direct 
detection from scattering experiments. 

There are several different types of detec- 
tor used in direct detection experiments. 
Fundamentally, WIMP detectors all work 
off the assumption that electrically neu- 
tral WIMPs will not interact with atomic 


electrons, but do have a very small chance 


of colliding elastically with atomic nu- 
clei. This collision would produce a de- 
tectable energy transfer. The summary be- 
low largely follows the information given 
in [8]. 

Anorganic Crystal Detectors were used in 
the very first experimental attempt at dark 
matter detection in 1987 [16]. This type 
of detector uses a high purity Germanium 
(HPGe) crystal as a target. For a semicon- 
ductor such as this, only a small amount 
of energy is required to form an electron- 
hole pair. A dark matter interaction should 
be easily detectable as many signal carri- 
ers will be produced in the event. The 
biggest issues with this detector are scal- 
ability and the lack of shielding to back- 
ground radiation. Germanium crystal be- 
comes very costly as the size of the target 
increases. This prevents making massive 
detectors of this kind which, given the rar- 
ity of dark matter interactions, significantly 
reduces the chances of a successful detec- 
tion. 

Cryogenic Detectors operate at tempera- 
tures of less than 50mK. They use dielectric 
crystalline targets such as Ge or Si. The 
system is weakly coupled to a heat bath 
and, upon interaction, two signals are mea- 
sured simultaneously: a phonon signal and 
an ionisation signal. However, operation at 
such low temperatures makes these detec- 
tors very expensive and limited in size. 


Bubble Chambers maintain superheated 


liquids at a temperature just below their 
boiling point. A dark matter interaction 
will cause a phase change in the liquid and 
a microscopic bubble of gas will form. A 
common problem with WIMP detection is 
the challenge of distinguishing nuclear re- 
coil from electronic recoil which both pro- 
duce similar signals. Electronic recoil can 
be caused by background y or 8 radiation. 
For a bubble to be formed, the energy de- 
posited within a given length (=) must 
exceed a critical value which varies with 
the temperature of the liquid. The amount 
of energy that an incident particle imparts 
onto the liquid depends on its stopping 
power. Thus, different particle interactions 
produce different ionisation densities which 
provides a means to discriminate between 
different particles. [17] Therefore, the ex- 
periment can be tuned so that bubbles can 
only be formed through interactions with 
neutrons, a@ particles or WIMPs. This elim- 
inates the possibility of eletronic recoil cre- 
ating bubbles. If the background radiation 
can be significantly reduced, neutrons and 
a particles would not be present in the de- 
tector. Therefore, the only remaining cause 
of bubble formation is WIMP interaction, 
meaning WIMPs can be unambiguously de- 
tected. One major benefit of bubble cham- 
bers is that they can achieve target masses 
of around 50kg, making them one of the 
leading prospects in dark matter detection. 


Noble Liquid Detectors use liquid Argon or 


Xenon as the target material. A diagram of 
this kind of detector is shown in Figure 3. 
Argon and Xenon are easily ionised. In the 
event of a WIMP collision, ions, X*, and 
excited atoms, X*, are produced. The ex- 
cited atoms combine with neutral atoms to 
form excimer states, Xš, which then decay, 
giving off scintillation light. An electric 
field, a “drift field”, is applied across the 
detector fluid to pull ions and ionised elec- 
trons apart to prevent them from recombin- 
ing. In order to distinguish electronic recoil 
from nuclear recoil, two signals (S1 and S2) 
are measured. Both of these signals are de- 
tected by photosensors, traditionally pho- 
tomultiplier tubes (PMTs). The ionised 
electrons provide the second signal. To pro- 
duce the S2 signal, a gas phase is positioned 
at the top of the detector. When the elec- 
trons arrive at the phase boundary, they are 
pulled across the gas phase by a strong “ex- 
traction field”. They then produce a signal 
proportional to the number of electrons ex- 
tracted. In electronic recoil, the number 
of ionised electrons is much greater than in 
nuclear recoil. Thus, a weaker S2 signal 
is expected for nuclear recoil. The ratio 


of S2/S1 can therefore distinguish between 


the two types of interaction. 


D. The DarkSide-20k Experiment 


One upcoming attempt at directly detect- 
ing WIMPs comes in the form of the 


Figure 3: The photomultiplier tubes (PMTs) first 
detect the scintillation signal (S1). The two electric 
fields then extract the electrons for the detection 
of the second signal (S2). [8] 


DarkSide-20k project. The experiment will 
run for a minimum of 5 years and aims to 
finally detect WIMP dark matter. It will 
use a Noble Liquid Detector with Liquid 
Argon (LAr) as the target material inside 
a Time Projection Chamber (TPC). 

As has been discussed, this type of de- 
tector can distinguish between nuclear and 
electronic recoil. Therefore, background y 
and 8 radiation is not of much concern. 
The most dangerous type of background 
radiation for this experiment is from neu- 
trons. Single neutron-nucleus interactions 
produce identical signals to WIMP-nucleus 
interactions. Shielding provides protection 
against neutrons entering the TPC from 
the external environment. However, the 
components within the TPC may emit ra- 


diogenic neutrons which cannot be shielded 


against. Neutron induced recoils such as 
this are suppressed by a liquid-scintillator 
neutron veto (LSV). This is a separate de- 
tector surrounding the TPC which detects 
internal and external neutrons. This allows 
recoils from these neutrons to be identi- 
fied and rejected. The LSV is surrounded 
by a Cherenkov detector called a Water 
Cherenkov Veto (WCV). This provides ad- 
ditional stopping power for external neu- 
trons and y radiation, as well as lower- 
ing the energy threshold for neutron detec- 
tion in the LSV. Furthermore, the materials 
used in the detector are carefully chosen to 
reduce radiogenic neutrons as much as pos- 


sible [18, 19]. 


E. Photomultiplier Tubes and 


Silicon Photomultipliers 


In the past, Noble Liquid Detectors have 
used photomultiplier tubes (PMTs) to de- 
tect the 2 signals given off by the target 
material [20, 21, 22, 23]. PMTs take ad- 
vantage of the photoelectric effect and elec- 
tron multiplication. A diagram of a PMT 
is shown in Figure 4. A photon strikes a 
surface with a low work function (a photo- 
cathode) and an electron is emitted. The 
electron is accelerated to a second surface 
where the impact causes the emission of 
several electrons. This process is repeated 
so that many electrons are available for sig- 


nal detection [24]. 
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Figure 4: A photon passes through a window and 
begins the process of electron multiplication. The 
electrons then provide a signal which indicates the 
incidence of the photon. [24] 


The DarkSide-20k project will use a sim- 
ilar, but more advanced technology: Sili- 
con Photomultipliers (SiPMs). These of- 
fer greater photon detection efficiency than 
PMTs, better single-photon resolution, and 
lower electrical power requirements. Fur- 
thermore, a key component of PMTs is 
glass which produces radiogenic neutrons. 
SiPMs contain no glass, making them a 
much better choice of photosensor for this 
project as they help to keep background ra- 
diation to a minimum. 

The following discussion of SiPMs uses in- 
formation given in [25]. The fundamen- 
tal technology at the heart of the SiPM 
is the avalanche photodiode (APD). The 
APDs are designed to work in Geiger- 
mode, meaning they are capable of detect- 
ing single-photons. An APD consists of a 
photodiode, a resistor and (although not 


Figure 5: Circuit diagram of a Geiger-mode APD. 
A photon (shown as hv) strikes the photodiode and 
current begins to flow down two paths: Ipes and 
Isap. [25] 


always) a capacitor as shown in Figure 5. 
The breakdown voltage of a photodiode is 
the point at which the reverse bias is broken 
and current flow rapidly increases. In order 
to operate in Geiger-mode, a bias voltage, 
Vias, Must be applied which is greater than 
the breakdown voltage, V;,. Initially, even 
with the application of Vias, there is no cur- 
rent flow in the photodiode. Upon the in- 
cidence of a photon, an electron-hole pair 
is created. The electron is accelerated by 
Voias and additional electrons undergo im- 
pact ionisation, causing avalanche multipli- 
cation. Now that the current is flowing, the 
voltage increases across the resistor. Across 
the photodiode the voltage now drops be- 
low V,,. This allows the directional bias 
of the photodiode to recover and stop the 
current flow. With no current, the voltage 
across the resistor drops back down again. 


Finally, the voltage across the photodiode 


rises back up to Vbias, making the system 
ready to detect the next photon. 

There are two key limitations of an individ- 
ual APD. Firstly, an APD is incapable of 
detecting how many photons are incident 
on the photodiode at one time. Secondly, 
the area of the photodiode is limited in size, 


for two reasons: 


e Thermal noise in the diode increases 
with active area (the area capable of 


detecting photons). 


e A larger active area increases the ca- 
pacitance of the APD, and a greater 
capacitance increases its recovery time 
(the time that the APD is insensitive 


to incoming photons). 


To work around these limitations, a SiPM 
is made up of many APDs (hundreds or 
thousands per square millimetre). In a 
SiPM, each APD is called a microcell. Fig- 
ure 6 shows an array of microcells forming 
a SiPM. 

A SiPM allows for the detection of multiple 
photons (as shown in Figure 7). The need 
for small-area photodiodes, as well as re- 
sistor and capacitor elements, means there 
is dead space in each microcell. This limits 
the total active area. The percentage of the 
total area that is active is called the fill fac- 
tor. A higher fill factor allows for greater 
photon detection efficiency (PDE) - a key 
parameter of SiPMs. 
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Figure 6: A SiPM made up of a 4x3 array of 
parallel-connected microcells. The bias voltage is 
applied across the cathode and the anode. [25] 


Other parameters are: 

Gain - the internal gain of the SiPM in re- 
sponse to a single-photon interaction. 
Dark Count Rate - the number of false sig- 
nals measured per second due to thermal 


noise. Dark count signals have the same 


amplitude as single-photon pulses. (See 
Figure 7.) 

Crosstalk - The output signal caused by one 
microcell randomly firing one or more ad- 
ditional microcells in the array. (See Fig- 
ure 7.) 

Afterpulsing - Secondary avalanche events 
that are caused by trapped charge carriers 
from a previous avalanche event. 

Dynamic Range - The ability of a SiPM 


to detect incoming photons without satu- 


ration. 
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Figure 7: An oscilloscope reading of a SiPM signal 
after low-level light pulse. The ability of the SiPM 
to detect multiple photons is shown by the 5 dis- 
crete voltage levels which represent the detection 
of 1 photon to 5 photons [25]. 


F. Crosstalk and Afterpulsing 
Probability 


Photon-detection and dark counts are 
known as primary events. When a micro- 
cell is triggered due to a primary event, 
it has the potential to generate secondary 
events (crosstalk and afterpulsing). 
Crosstalk occurs from a primary event 
when the primary microcell simultaneously 
triggers one or more of those immedi- 
ately next to it. Microcells triggered via 
crosstalk are capable of producing fur- 
ther crosstalk events in additional micro- 
cells. Thus, crosstalk is best described as 
a branching Poisson process, shown in Fig- 
ure 8. [26] 

For such a process, the probability, P(x), of 
x events occurring is given by a generalised 


Poisson distribution with a shape charac- 


Random primary (Photo) events 
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Figure 8: An initial photo event triggering a micro- 


cell which then causes a chained process of branch- 
ing crosstalk in neighbouring microcells [26]. 
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terised by the parameters 0 and A [27]: 
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Figure 9 shows the shape of the generalised 
Poisson distribution, plotted using pyplot, 
for 0 = 3.5 and A = 0.51. The choice for 
0 here is fairly arbitrary, but A has been 


measured as ~ 0.51 and so was set to this 


value. [26] 
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Figure 9: The shape of the generalised Poisson dis- 
tribution with 0 = 3.5 and à = 0.51. The probabil- 
ity of observing a given number of crosstalk events 
is plotted against the number of crosstalk events. 


When a primary event occurs, the cur- 


rent in the microcell surges before the de- 


vice self-regulates to halt the current flow. 
During this process, some charge carriers 
may become trapped in the device before 
discharging at a later time to produce a 
second pulse. This kind of event is an af- 
terpulse and has an amplitude less than or 
equal to that of a primary pulse. [28] The 
probability of an afterpulse occurring at a 
given time is described by a Poisson distri- 
bution. [29] At a time, At, after a primary 
pulse, the probability, P,,(At) of an after- 
pulse taking place is given by: 


Pap(At) = e7% (2) 
Tap 

Where Tap is the time constant. The after- 
pulsing process could be better described 
using two time constants, one short and 
one long. However, since there is usually 
a time window applied for the detection 
of secondary events, those corresponding to 
the longer time constant are not detected. 
Thus, only one time constant is needed. 
The value of Tap can range from 10ns to 
10us. [29] For the Hamamatsu SiPM, Tap 
has been measured to be ~ 15ns. [30] The 
shape of the afterpulsing probability distri- 


bution for Tap = 15ns is shown in Figure 10. 


When analysing the signal from a SiPM, 
the different noise types can be distin- 
guished by plotting the amplitude of each 
pulse against the time since the one that 


preceded it. [28] An example of this is 
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Figure 10: The probability of an afterpulse occur- 
ring is plotted against the time since the original 
primary event with Tap = 15ns. 
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Figure 11: The different SiPM noise types appear 
in different areas on a plot of amplitude vs time 
since preceding pulse. [28] 


shown in Figure 11. There are several im- 
portant features visible in this plot. Firstly, 
dark counts make up a thick horizontal line 
with an amplitude of 1.0. 

Afterpulses are visible with amplitudes in 
the range of 0.5 to 1.0 occurring at times of 
50 to 200ns. The reason for this amplitude 
range is that a minimum amplitude of 0.5 
was required for detection. 

Finally, direct crosstalk can be seen above 
the dark count line, and delayed crosstalk 
can be seen to the left of it. 


This plot from [28] demonstrates how the 
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Figure 12: The model of a SiPM circuit used as a 
basis for the first stage of this project. [28] 


different noise types in a SiPM can be dis- 
tinguished from each other, thus enabling 


characterisation of the device. 


II. Methods 


A. Circuit Design 


The first stage of this project was to build 
a model of a SiPM in circuit simulation 
software. The software used was Multisim. 
The fundamental circuit design (shown in 
Figure 12) was based off of the one given 
in [28]. The design includes a SiPM (an ar- 
ray of microcells as in Figure 6) connected 
to an operational amplifier (op-amp). An 
oscilloscope is used to read the output of 
the circuit. There are two reasons for the 
use of the op-amp. The first is to match 
the impedances of the SiPM circuit and 
the oscilloscope to prevent signal loss. Sec- 
ondly, the op-amp translates the current 
signal from the SiPM to voltage which can 
be read on the oscilloscope. From this ba- 


sis, the model was further developed to 


match the specifications of three different 
SiPMs: SensL C-Series, Hamamatsu, and 
FBK NUV. The specifications for the SensL 
C-Series and Hamamatsu SiPMs are given 
in [31]. Those for the FBK NUV SiPM were 
given in [18]. 


B. Modelling a SiPM Pulse 


Despite this project being in three distinct 
stages, only the first was expected to be 
completed. The second stage of this project 
was developing a program which could real- 
istically simulate a signal from a SiPM with 
all the additional noise types (dark counts, 
crosstalk, and afterpulsing). The first step 
was to input a pulse shape which replicated 
that of a SiPM. Mathematically, this shape 
can be described by a function with three 
time constants Tr, Ts, and 7). [32, 33] (Here, 
r, s, and l stand for rise, short, and long re- 
spectively.) The voltage output, V, of the 
SiPM as a function of time, t, is shown be- 


low: 


Where to is the onset (i.e. the start time of 
the pulse) and a is a scaling factor which 
allows for control of the magnitude of the 
pulse. Using this equation, a primary pulse 
representing a photo event can be created 
which always has the same onset and am- 
plitude. Noise can then occur from (or ran- 


domly after) this pulse. When program- 
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ming the pulse function, tọ was set to 50ns 
and a was set to 10. The three time con- 
stants, Tr, Ts, and 7 were set to Ins, 6ns, 


and 150ns respectively 


C. Modelling Dark Count 


The first noise type added to the signal was 
dark count. The number of dark count 
events in a given time window follows a 
Poisson distribution. [31] The mean value 
of this distribution is given by the product 
of the size of the time window and dark 
count rate (DCR). Thus, the amount of 
dark count can be adjusted via a DCR pa- 
rameter. This parameter is set to a typical 
value 0.002 which corresponds to a DCR of 
2MHz. 

In the simulation, dark counts occur in- 
dependently of the original pulse and the 
number observed in a given time-frame is 


chosen from the Poisson distribution. 


D. Modelling Crosstalk 


The probability distribution given in Equa- 
tion (1) was used to simulate crosstalk 
pulses. The value of A for crosstalk has 
been measured to be ~ 0.51. [26] This dis- 
tribution was programmed such that the 
two parameters, 0 and A could be altered. 
The value of A was set to 0.51. An ap- 
propriate initial value of 6 was deemed to 
be the value that would give the probabil- 


ity of having no crosstalk as 0.5. Using 
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Figure 13: The probability distribution of crosstalk 
with 0 = In(2) ~ 0.693 and à = 0.51. 


Equation (1) and setting à = 0.51, z = 0, 
and P(x) = 0.5 it can be calculated that 
6 = In(2) ~ 0.693. The output of this func- 
tion is shown in Figure 13. The next step 
was to program a function which would ex- 
tract a random number from this distribu- 
tion. To verify the accuracy of this func- 
tion, a million random choices were carried 
out. The shape of the distribution was re- 
produced by these choices which confirmed 
the function was accurate (see Figure 14). 
This was then incorporated into the pro- 
gram. 

When running through a simulation, at 
each primary event, a random number of 
crosstalk events is chosen to occur from the 
generalised Poisson distribution. The time 
at which these occur after the preceding 
primary pulse is then chosen from a Pois- 
son distribution with a mean of 50ns (an 
appropriate value for crosstalk [28, 34]). A 
pulse with amplitude proportional to the 


number of crosstalk events is produced in 
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Figure 14: The shapes of both the calculated prob- 
ability distribution and the distribution built up 
from a million random choices. 


the signal output. The maximum number 


of crosstalk events is set at 5. 


E. Modelling Afterpulsing 


Afterpulsing requires the trapping of 
charge carriers from a primary pulse. The 
number of carriers trapped follows a Pois- 
son distribution, [30], and the time at which 
they discharge also follows a Poisson distri- 
bution as given in Equation 2 [29]. The 
simulation of afterpulsing occurs in two 
steps. Firstly, for each primary (photo or 
dark) event, a number of trapped charge 
carriers is chosen from a Poisson distribu- 
tion with a mean of 1. If this number is 
zero, no afterpulse occurs. If it is non-zero, 
then the afterpulse occurs at a time chosen 
from the Poisson distribution in Equation 2 
(with a time constant of Tap = 500ns). 


In reality, the amplitude of an afterpulse 


is less than or equal to that of a primary 
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Figure 15: The circuit model of the Hamamatsu SiPM. An array of microcells is connected to an op-amp 
with the signal being read through an oscilloscope. 


Circuit Quench Resistance (kQ) | Microcell Capacitance (fF) | Feedback Resistance (kQ) 
Hamamatsu 211 3 200 
SensL C-Series 633 3 200 
FBK NUV 2000 40 1 


Table 1: The strengths of the microcell resistor (quench resistor), the microcell capacitance and the 
feedback resistance are displayed for each model circuit. [18, 31] 


pulse. [28] However, for this simulation, the III. Results 

afterpulsing amplitude was chosen to be 

half a primary pulse so that they could be A. Building a Model SiPM Circuit 

easily distinguished from dark counts. 

E E a E E Three different models of SiPMs were built. 
These models were based on the SensL 


C-Series, Hamamatsu, and FBK SiPMs. 


SiPM output results in a complex signal 


that can be used to develop software which 


Each ] has th truct dif- 
can accurately and efficiently analyse such D BOS ANG a 


fering only in the resistances and capaci- 
pulses. 


tances of the components. The model of 
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9.157V 


| 


m 
152ns 


Figure 16: The output of the Hamamatsu SiPM 
circuit. The pulse has a peak voltage of 9.157V 
and it takes 152ns to recover to its original value. 


the Hamamatsu SiPM is shown in Fig- 
ure 15 and its output in Figure 16. Each 
circuit has four distinct sections. Firstly, 
there is the pulse generator (left-hand side 
of Figure 15, labelled ‘V3’) which produces 
a square-wave of 1V amplitude with a pe- 
riod of lus (replicating a photon detection 
every microsecond). This signal is then fed 
into an array of 8 microcells. This consti- 
tutes the SiPM part of the circuit which re- 
stricts the current flow, producing a pulse 
signal. The third section of this circuit is 
a transimpedance amplifier with a 200kQ 
feedback resistor. The amplifier translates 
the incoming current into voltage so that 
it can be read on an oscilloscope (the fi- 
nal component). The different properties 
of each circuit model are given in Table 1. 
These models have not been built simply to 
verify that this circuit accurately produces 
realistic SiPM pulses. They are effectively 
a toolkit with which the factors affecting 


the properties of a SiPM can be investi- 
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Slope = (1.01 + 0.03)x10°? so! * 
A 


P 


0 200 400 600 


Microcell Resistance (kO) 


800 


Figure 17: Microcell resistance is plotted against 
the value of 7 and a linear relationship is found. 


gated. The components of the circuit can 
be altered as desired to probe the system’s 
behaviour. This is demonstrated by alter- 
ing the values of all the microcell resistors 
in the Hamamatsu model and observing the 
effect this had on the time constant 7;. The 
oscilloscope data from Multisim was input 
into a program where a curve-fitting func- 
tion fit the data to the model in Equa- 
tion 3. This function then output the fit 
parameters such as 7;. The resistances were 
changed from 100kQ to 900kQ. Figure 17 
shows the linear relationship between these 
two parameters. 


B. Simulating a SiPM Signal 


Using the models for SiPM noise previously 
given, a realistic simulation of a SiPM was 
developed. Three noise types (in addition 
to some thermal noise) were added to the 
pure output from the circuit. These were 


dark count, crosstalk, and afterpulsing. 


1000 


SiPM pulse 
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0 250 500 750 1000 1250 1500 1750 2000 
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Figure 18: The output of the simulation with only 
thermal noise applied to the original photo pulse. 
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Figure 19: A signal with the original photo pulse 
as well as two dark pulses, one at 1068ns and the 
other at 1542ns. 


The shape of the signal with no complex 
noise is shown in Figure 18. The model for 
dark count was added and an example of 
the resultant signal is given in Figure 19. 
Next, the crosstalk was added which pro- 
duced signals such as that in Figure 20. 
The final simulation of the signal with dark 
counts, crosstalk, and afterpulsing is shown 


in Figure 21. 


15 
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Figure 20: Two double-dark pulses occur at 948ns 


and 1699ns. After each of these there is a crosstalk 
pulse, these are at 1207ns and 1899ns. 
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Figure 21: The full simulation of a signal from a 
SiPM. This includes each type of noise that can 
occur. 


C. Signal Analysis 


In order to analyse the simulated signals, 
a peak finder function was employed. This 
function detects all the peaks in a signal 
and returns the amplitude of each peak and 
its distance from the one before it. To find 
the best parameters for the function, a loop 
was programmed to run through over a mil- 
lion configurations. The optimal parame- 


ters were found and incorporated. The re- 


sult of applying this function to a signal 


Filtered signal with peak finder 
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8 8 8 & 8 
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Figure 22: The result of running a peak detection 
program over the SiPM signal in Figure 21. It is 
clear that the peaks are effectively detected. 


Relative ampliude 


Time distance (Nanoseconds) 


Figure 23: The amplitude of every peak from 1343 
signals is plotted against the time since the peak 
that preceded it. 


window is shown in Figure 22. 

After building the peak finder, many sig- 
nals could then be analysed in one go. This 
made it possible to build up a large amount 
of information about the behaviour of the 
SiPM model. The amplitude of each peak 
was plotted against its distance in time 
from the peak before it. The result is shown 
in Figure 23. The key features of this plot 


are: 


e The horizontal line at 0.5 amplitude 


made up of afterpulses. 
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e The horizontal line at 1.0 amplitude 


made up of dark counts. 


A vertical bar between 40 and 60 
nanoseconds which is made up of 
crosstalk pulses occurring at their 
most probable time. The number of 
pulses decreases as the amplitude in- 
creases since larger numbers of simul- 


taneous events are less likely (see Fig- 


ure 13). 


The diffuse scattering of points above 
the dark count line due to crosstalk 


events at less probable times. 


The thin line at ~ 0 amplitude due to 


anomalous detections of thermal noise. 


The simulation was programmed such that 
the different noise types could be con- 
trolled. For example, the amount of 
crosstalk can be altered by changing the 
value of 0 in Equation 1, or dark count can 
be varied by changing the DCR parameter. 
The amount of crosstalk was significantly 
reduced by setting 0 = 0.2 and the analysis 
program was used to make a second scatter 
plot. This is displayed in Figure 24 which 
clearly shows that the amount of crosstalk 
is reduced. Furthermore, the amount of af- 
terpulsing and dark count on the plot has 
not changed. 

Finally, a third plot (Figure 25) was made 
with a greatly increased crosstalk probabil- 


ity, but a much smaller DCR. The value of 0 


Relative ampliude 


Time distance (Nanoseconds) 
Figure 24: 1343 signals were analysed to produce 
a second plot. The value of 0 was changed to 0.200 
and the resultant reduction in crosstalk is clearly 
visible. 


Relative ampliude 


Time distance (Nanoseconds) 
Figure 25: For the final plot, the value of 0 was 
raised to 2 and the DCR parameter was lowered to 
0.0005. The reduction in dark count and increase 
in crosstalk are visible. 


for this data-set is 2, the maximum number 
of crosstalk events was increased to 10, and 
the DCR parameter was set to 0.0005 which 
corresponds to a DCR of 500kHz. Figure 25 
shows the reduced dark count (and there- 
fore also reduced afterpulsing) as well as 


much more crosstalk. 


IV. Discussion 


It is important to note that sections B, C, 


D, and E in Methods are actually part of 
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the results of this project. Even though 
they are part of what was achieved, they 
have been included in the Methods section 
because they more appropriately describe 
how. In other words, the modelling and 
simulation of each noise type were results 
on their own, but they better describe how 
the SiPM simulation was built. 

The results of this project have many po- 
tential applications both in the DarkSide- 
20k project and in wider researh of SiPM 


technology. 


A. The Model SiPM Circuits 


The original goal for this project was sim- 
ply to build a circuit that could replicate 
the behaviour of a SiPM. This circuit was 
created, but more importantly, three dif- 
ferent SiPM models were developed from 
this basis. These were models of the SensL 
C-Series, Hamamatsu, and the FBK NUV 
SiPMs. They differed only by the strengths 
of the resistors and capacitors in the cir- 
cuits. Through alteration of these proper- 
ties, the relationship between quench resis- 
tance and 7; was found. The models with 
lower quench resistance have a shorter re- 
covery time. This demonstrates how these 
models make up a toolkit for further inves- 
tigation into how circuit properties affect 
the behaviour of a SiPM. Thus, what has 
been produced in this project has the po- 


tential to aid the development of more effi- 


cient SiPMs in the future. In spite of this, 
the curve-fitting function used to evaluate 
Tı was not able to reliably calculate values 
for T, and Ts. The rise times of the pulses 
can be in the order of picoseconds. [34] 
Therefore, it is likely that the surge of cur- 
rent during a primary event happens too 
quickly for the program to accurately de- 
tect. This means the data which encodes 
the information about 7, is not precisely 
recorded, preventing 7, from being accu- 


rately calculated. 


B. SiPM Output Simulation 


A simulation was created of a SiPM signal 
output with the original photo pulse and 
all the noise that can stem from or accom- 
pany it. This simulation is not only a real- 
istic representation of a SiPM signal but it 
is also controllable. The levels of each type 
of noise can be raised or lowered by altering 
the parameters of the governing equations. 
Despite being a reasonable output from 
a SiPM, the simulation isn’t without its 
flaws. Firstly, the afterpulsing is simplified. 
It has been previously shown that the am- 
plitude of an afterpulse increases with the 
amount of time since the pulse that pro- 
duced it (see Figure 11). [28, 34] In this 
simulation however, all afterpulses have an 
amplitude of 0.5. Furthermore, the distri- 
bution of crosstalk in time does not agree 


with that in [28] or [34]. 
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As has been demonstrated in this paper, 
this simulation provides the basis from 
which software can be designed to analyse 
the simulated output and, by extension, the 


output from any real SiPM. 


C. Signal Analysis 


The signal analysis software relies on a peak 
finder function. It is possible that this func- 
tion will not resolve two or more separate 
peaks if they are too close together. The 
efficiency of this function is calculated as 
the number of peaks detected divided by 
the total number of events that occurred. 
The value was estimated by collecting data 
from 200 SiPM signal simulations. When 
multiple crosstalk events occur, they hap- 
pen simultaneously. This means multiple 
events are read by the peakfinder as only 
one event. When analysing SiPM signals, 
the number of crosstalk events is revealed 
by the amplitude of the peak. This makes 
it acceptable to merge multiple crosstalk 
events into one peak. However, when cal- 
culating the estimate, this results in many 
events being ‘missed’, even though this is 
not actually the case. Whilst there is a 
small number of ‘extra’ detections due to 
thermal noise, these are outweighed by the 
number missed due to the simultaneity of 
crosstalk. Therefore, this estimate provides 
a lower limit on the function’s efficiency. It 


was found that the peak finder has an ef- 


ficiency of at least 80.23%. The true effi- 
ciency is higher than this value, but it needs 
to be made more efficient for testing SiPMs. 
The whole purpose of this project was to 
prepare for the arrival of the DarkSide-20k 
SiPMs at Warwick, where they will un- 
dergo room-temperature testing. The sig- 
nal simulation provided the foundation to 
build the analysis software to be used for 
this testing. This means that when the 
SiP Ms arrive, their behaviour can be inves- 
tigated and their properties characterised. 
The levels of noise from all the SiPMs can 
be measured so that once they go into the 
detector, an expected amount of noise will 
be known and can therefore be disregarded. 
These tools do not only apply to DarkSide- 
20k, but also apply to any attempts at 


SiPM characterisation. 


V. Conclusions 


The work performed throughout this 


project was done in three stages: 


e Model circuits were created. 


e A complex output from a SiPM was 


simulated. 


e Software was written to analyse SiPM 


signals. 


Model circuits were built which accurately 
and reliably reproduce the behaviour of ex- 
isting SiPMs. These can be used for fur- 
ther investigation of SiPM properties and 
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efficiency. 

Three types of complex noise were simu- 
lated: dark count, crosstalk, and afterpuls- 
ing. This produced a realistic SiPM output 
signal which could then be analysed. How- 
ever, there were issues with this simulation 
which reduced its ability to provide a com- 
pletely accurate representation of a SiPM. 
A higher level of realism will aid the im- 
provement of the analysis software as any 
issues that will occur with real SiPMs, will 
also occur with the simulation. 

Analysis of this signal was performed which 
demonstrated how the different noise types 
in a SiPM can be identified during the test- 
ing process. 

There are numerous opportunities for fur- 
ther research from this project. With more 
time, the curve-fitting function would be 
improved and made more reliable. The 
model circuits would also be used to ex- 
tensively investigate how each component 
of the SiPM circuit affects its overall be- 
haviour. Also, the simulation of the after- 
pulsing noise would be developed to repli- 
cate the behaviour of real SiPMs such as 
those studied in [28] and [34]. The ideal 
parameters of the afterpulsing probability 
equation would be found such that they 
produce pulses within 50 to 200ns of the 
preceding pulse. This would be a much 
more realistic time-distribution of after- 
pulsing. An extra level of complexity with 


afterpulsing, which was not included in this 


simulation, is the fact that the pulse am- 
plitude scales with the time since the pulse 
that generated it. Further development of 
this work would add this complexity. If ad- 
ditional time was available, the efficiency of 
the peak finder function will be more pre- 
cisely calculated. Only a lower limit could 
be applied in this project however, after 
more work, it is possible a more precise 
method could be developed. Finally, the 
analysis software developed in this project 
will be used to test and characterise the 
SiPMs for use in the DarkSide-20 collabo- 


ration. 
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